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ABSTRACT 

Expanding blast waves are ubiquitous in many astronomical sources, such as supernovae remnants 
(SNRs), X-ray emitting binaries (XRBs) and gamma-ray bursts (GRBs). I consider here the dynamics 
of such an expanding blast wave, both in the adiabatic and the radiative regimes. As the blast wave 
collects material from the surrounding, it decelerates. A full description of the temporal evolution of 
the blast wave requires consideration of both the energy density and the pressure of the shocked ma- 
terial. The obtained equation is different than earlier works in which only the energy was considered. 
The solution converges to the familiar results in both the ultra-relativistic and the sub-relativistic 
(Newtonian) regimes. 

Subject headings: hydrodynamics — ISM: jets and outflows — relativity — shock waves 



1. INTRODUCTION 

Expanding blast waves are one of the most common 
phenomena in many astronomical transients. Inter- 
action of the expanding shells of supernova remnants 
(SNRs) with their e nvironment was studie d more than 
three decades ago (| Chevalier! 119761 I1982T) . Similarly, 
in X-ray emitting binaries (XRBs), mildly relativistic 
(Lorentz factor F ~ few) expandi ng radio blobs are ob- 
serve d for nearly two decad e s now (Mirabel & R odriguez 
[1991 IHiellming k Rupenl IJ995t iFender et al.l I1999JT 
In gamma-ray bursts (GRBs) relativistic blast waves 
are an inherent part of the GRB "fireball" model 
(IGoodmanlfrOSl IPaczvnskil Il98l IShemi k Piranl Il990t 
iRees k Meszaroa [l992: Naravan et al.lfl992h . and are the 
source of the afterglow emission frequently seen. 

These expanding blast waves originate from a stel- 
lar explosion (such as in GRBs or supernovae), or a 
rapid ejection of material (as in XRBs). As they prop- 
agate through the ambient medium, they collect mate- 
rial and decelerate. The expansion may be adiabatic, or 
highly radiative. The dynamics of the blast wave expan- 
sion in b oth these scenarios were extensively studied in 



1997; 



the past [ Blandford k McKee 1976 ;lKatz k Pira: 
Chia ng k Dermerl ll999; Piran 1999; [Huang et al 
van Paradijs et aLl i2000)rl These formulations are used 



in 
il. 



1999; 
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McKec ( 1976) provided a full solution in 
the ultra-relativistic, adiabatic limit, the original equation that 
describes the blast wav e evolution in t he rad iative scenario was 
written and solved by Katz & Piran| 1119971 ). It was general- 
ized to include b o th th e radiat i ve an d adiabatic scenarios by 
Chiang & Dcrmcr (1999): Piran (1999). However, it was later 
shown by IHuane et al.lill999l) that this equation does not con- 
verge to the correct asy mptotic solution for adiabatic expansion 
in the Newtonian limit. Huang et al. (1999) thus modified this 
equation, solved the revised equation and showed that its solution 
converges to the correct asymptotic solutions in both the ultra- 
relativistic and the Newtonian limits, for both ad iabatic and ra- 
diative scenarios. Bianco & Ruffini (2004, 2005a . b), apparently 
being unaware of the work of Huang et al. (1999), have re-solved 
the dynamical equation written by Chiang & Dcrmcr ( 1999); Piran 
(1999), both numerically and analytically. Apart from the work by 
Blandford & McKec (1976), all the other works relied on similar 
basic assumptions, in which the contribution of the shocked-ISM 
pressure was neglected, and are therefore conceptually incorrect. 



until today as a basis for calculating the expected radia- 
tion during the deceleration phase of the blast waye evo- 
lution (e.g.. iNaravan et all 12011b iGranot k Pi7anll2012t 
Ivan Eerten k MacFadyeiJl2012t iShen k Matznerll2012h . 
The dynamical calculations in the above mentioned 
works are based on a "toy model" , in which a basic as- 
sumption is that the interaction between the blast wave 
and the interstellar material (ISM ) can be described 
by a series of ine lastic collisions (jKatz k Piranl 11997b 
iHuang et al.l[l9 99). While this prescription asymptotes 
to the known analytical solutions in th e ultra-relativistic 
limit L > 1 dMeszaros k Reeslfl997J ) and in the non- 
relativistic limit (|Sedovlll959D . it neglects the contribu- 
tion of the swept-up material to the internal pressure. 
This, in turn, affects the evolution of the blast wave dy- 
namics. In this letter, I revise the basic assumptions of 
the blast wave - ISM interaction scenario, and re-derive 
the equations that govern the blast wave evolution. As 
will be shown below, the obtained equations are different 
than the ones previously used. Nonetheless, they asymp- 
tote to the known results at the ultra-relativistic as well 
as the sub-relativistic limits. 

2. DYNAMICAL MODEL 

Consider an explosion that ejects mass M and creates 
a shell that propagates into the cold ISM. In front of 
the expanding shell is a blast wave. The system is thus 
composed of 3 regions: (1) the unshocked ISM, (2) the 
shocked ISM, and (3) the ejected shell material. For 
simplicity, I assume that a reverse shock wave (if cre- 
ated) had long passed, hence the ejected shell material 
is cold. Furthermore, I assume that the thermodynami- 
cal quantities of the gas: rii,pi and ei (particle number 
density, pressure and internal energy density) are steady 
in each region. The pressure in each region is given by 
p' = (7 — l)(e' — p'), where p' = n'm p c 2 is the rest mass 
density and 7 is the adiabatic index. Here and below, 
primed quantities are in the comoving frame, and un- 
primed quantities are in the observer's frame in which 
the unshocked ISM is at rest. 

In the following , I will compare the results derived here to the re- 
sults derived by Huans et al. (1992), that were shown to converge 
to the correct asymptotic limits. 



The evolution of the blast wave as it propagates 
through the ISM is calculated by considering energy con- 
servation in the observer's frame. Let us assume that at 
time t the plasma propagates with Lorentz factor Fo Ne- 
glecting radiative losses and assuming that the gas can 
be described as a prefect fluid, the energy density in re- 
gion (2) as is measured by a distant observer is given 
by 



e 2 = T 00 = (e 2 + p' 2 )u°u () + p' 2 g m 

= [ 7 r 2 - (7 - 1)] 4 + ( 7 - i)(i - r 2 Km p c 2 

(1) 

Here, T 00 is the 00 component of the stress-energy ten- 
sor. Due to Lorentz contraction, the comoving volume of 
region (2) is V — TV. The total energy of the gas con- 
tained in this region (in the adiabatic case), as is viewed 
by a distant observer is thus 



eW 



( 7 - l)] 5 ^- + (7 - 1)(1 - T 2 ) n ^ cV 
( 7 -l)(l + r/3 2 )]7V 2 m p c 2 



E 2 (ad.) = [ 7 L 2 

= [7r 2 

where N% — n 2 V is the number of particles in region 
(2), f3 = (1 — T -2 ) 1 / 2 is the normalized bulk velocity 
of the plasma in this region, and use was made in the 
relation e 2 /n 2 = rm p c 2 , which is exact for any value 
of r as long as the unshoc ked ISM in region (1) is cold 
(|Blandford fc McKedll976ft . 

The calculation in equation [2] assumes no radiative 
losses. In order to allow the possibility that part of the 
thermal energy gained by the ISM as it crosses the shock 
wave is radiated, equation [2] is modified as follows. The 
energy calculated in equation [2] is the sum of three sepa- 
rate components: (1) the rest mass energy of the shocked 
ISM; (2) its kinetic energy; and (3) its thermal energy. 
The first two components sum up to riV 2 TOpC 2 . Only the 
energy in the third component can in principle be radi- 
ated. If a fraction e of the thermal energy is radiated, 
then the energy of the gas in region (2) is given by 



E 2 = {r + (1 - e)[jT 2 - r - (7 - 1)(1 + T(3 2 )}} toc 2 , 

(3) 
where to = -/V 2 to p is the mass of the shocked ISM. The 
material in region (3) is assumed to be cold, and so its 
(kinetic + rest mass) energy is E% = TMc 2 . 

A differential equation for the evolution of the plasma 
velocity is derived in the following way. Between times t 
and t + 6t, the plasma propagates a distance /3c6t, and 
an ISM of mass dm crosses the forward shock and gains 
kinetic and thermal energy. A fraction e of the gained 
thermal energy is assumed to be radiated, hence the radi- 
ated energy is 5E rad = e[7L 2 -L- (7- I)(l + L/3 2 )]dmc 2 . 
As it collects material, the plasma decelerates; at time 
t + St its Lorentz factor is L — dT (corresponding velocity 
(3 — d(3). Conservation of energy at times t and t + St is 



3 Note that the calculation here is general, and is not limited to 
ultra-relativistic speeds, where r>l. 



written as 

{r + (1 - e)[7r 2 - r - (7 - i)(i + rp 2 )}} m + tm 

= -dm + {(T - dT) + (1 - e)[ 7 (L - dT) 2 - (L - dT) 
-(7 - 1)(1 + (L - dT)(0 - dp) 2 )]} (to + dm) 
+(T - dT)M + e[ 7 (r - dT) 2 - (L - dT)- 

(7 - 1)(1 + (T - dT)(P - d/3) 2 )]dm. 

. W 

Re-arranging the terms in equation [4j it can be written 

as a dynamical equation for the evolution of the bulk 

motion Lorentz factor, 



dT 

dm 



7 (r 2 - 1) - (7 - i)r/3 2 



M + em + (1 - 6)to[2 7 L - (7 - 1)(1 + L" 2 )] 



( 5 ) 
Equation [S] thus describes the evolution of the plasma 

Lorentz factor due to its interaction with the ISM. This 

equation holds both in the ultra-relativistic (r> 1) as 

well as the sub-relativistic (/3 <C 1) lim i ts. It can be 

compared to equation (7) in lHuang et al.l (|1999ft . which, 

as described above, was derived based on the assumption 

of continuous inelastic collisions between the blast wave 

and the ISM, and hence does not contain the contribution 

of the shock-heated ISM to the pressure in region (2)0 

2.1. Asymptotic limits 

It is possible to obtain analytical solutions to the dy- 
namical equation [5] in the limits of ultra-relativistic and 
sub-relativistic limits. These are useful for comparison 
with former results. 

In the adiabatic scenario, e = 0, and equation[5]reduces 
to 

d£_ = 7 (r 2 - 1) - (7 - l)L/3 2 

dm 



M 



(6) 



to[2 7 L- (7- 1)(1 + L- 2 )] 

In the ultra-relativistic limit, r> 1, equation [5] can be 
re- written as 

d£_ 

dm 



(r 2 - 1) 



(M/ 7 ) + 2Lto' 



(7) 



Denoting by Tq the initial Lorentz factor of the flow, the 
evolution of the plasma Lorentz factor can be separated 
into two regimes. (1) Initially, M/j ^> Loto. In this 
regime, L ~ Tq. (2) For Tm 3> M/j one obtains the 
familiar solution, F oc to -1 ' 2 , which, for constant den- 
sity ISM (n oc r , to oc r 3 ) results in the well known 
solution F oc r~ 3 / 2 . I nterestingly, e quation [7] is similar 
to equation 8 of Huan get al.1 ([1999), with M replaced 
by M/7 ~ (3/4)M. This discrepancy has only a minor 
effect on the blast wave evolution. 

On the other extreme, the sub-relativistic limit /3 <C 1, 
equation |6] reduces to 



dm 

dm 



P 



M 



i[2-py 



(8) 



which admits the solution (3 oc m" 1 / 2 , for to ^> M. 
Thus, the general solution describing the blast wave evo- 
lution in the adiabatic scenario (equation [6]) for Tto 3> 
M/7 can be written as Tj3 oc m~ x l 2 . 

4 Note that equation (7) of IHuang et al.1 1119991 ) is obtained by 
setting 7 = 1 in equation [5] 

5 Equation [7] does not admit a third regime, Y^m >M> Tm. 



In the radiative scenario, e = 1. In the ultra-relativistic 
limit r>l, the dynamical equation (eq. [5]) becomes 



dT ^ jT 2 

dm M + to 



(9) 



Initially, M 3> to, resulting in a steady Lorentz factor, 
r ~ Tq. At a later stage, Tqitl 3> M and one obtains 
the decay law T(m) ~ M/jm (as long as to <C M). For 
constant density ISM, n oc r°, this leads to the familiar 
decay la w, T oc r~ 3 . This result is different than the 
result of iHuang et~aH (1999) by a factor (7)" 1 ~ 3/4. 
In the sub-relativistic limit, ^<1, equation[5]becomes 



<*(r/3) 

dm 



P 



M + m 



(10) 



with the solution /3 oc m _1 (for to 3> M). This so- 
lution is similar to the classical "snowplow" evolution of 
an expanding supernova remnants in the radiative regime 
IISmtzerlll968h . 

3. NUMERICAL SOLUTION 

A full solution of the dynamic equation [5] can easily 
be obtained numerically. In solving this equation, one 
first needs to determine the value of the adiabatic in- 
dex 7, which depends on the gas temperature. In cal- 
culating 7, I assume that the gas in region (2) main- 
tains a Maxwellian distribution with normalized temper- 
ature = \zbT jm v c 2 . The average energy per parti- 
cle in this region is thus given by < e' 2 /n' 2 m p c 2 >— 
K 1 (6- r )/K 2 (e- 1 ) + 3e. Here, ifr , K 2 are modified Bessel 
^-functions of the second kind (|Lightman et al.l 119751 
5.34). 

The ratio < e 2 /n' 2 m p c 2 >= T is determined by the 
shock jump conditions, and is thus known at any given 
instance. For a given T, a good fit to the normalized 
temperature is 



3 



r/3 + 1.07 (T/3) 2 

i + rp + i.07{r[3y 



(ii) 



This fit asymptotes to the exact solution in the limits 
r » 1 and )3 < 1. The maximum error found is less 
than 3 x 10~ 3 , for T/3 ~ 1. Once the gas tem perature 
is calc ulated, I use the polynomial fit given by iServicd 
(dllg ), to calculate 7 ~ (5 - 1.21937.? + 0.18203z 2 - 
0.96583z 3 + 2.325132 4 - 2.39332z 5 + 1.07136z 6 )/3, where 
z = 6/(0.24: + 9). This fit is accurate to 1CT 3 . 

Equation [5] is solved using 4th order Runge-Kutta 
method. I consider two sets of parameters, one repre- 
senting GRB and one XRB. In both scenarios, the blast 
wave is assumed to propagate into a constant density 
ISM, and hence the collected mass is 



dm — Aur nm p dR, 



(12) 



where n is the number density of the ISM and m p is 
the proton rest mass. Photons emitted as the plasma 
propagates a distance dR are observed at time dt, given 
by 

dR = r/3c (r + r/3) dt. (13) 

While equation [T3] is derived under the assumption that 
the observed photons are emitted from a plasma which 
propagates towards the observer, a more comprehensive 
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Fig. 1. — Temporal evolution of the bulk Lorentz factor, r. Thick 
lines are for "GRB" scenario (T = 278, M = 4 X 10 28 gr), while 
thin lines are for "XRB" scenario (To = 3, M = 3 X 10 23 gr). In all 
cases, ISM density n = 1 cm -3 assumed. Solid: adiabatic, t = 
(blue, black: GRB, XRB respectively); Dashed: radiative, e = 1 
(red, magenta: GRB, XRB respectively). 

calculation which considers the integrated emission from 
different angles to the line of sight results in a similar 
solut i on, up to a numeric al factor of a few (jWaxmanl 
IT997T: iFVeT fc Wii er d l2006h . 

The evolution of the Lorentz factor T, the momentum 
T/3 and the radius R are presented in figures Q] - [3] In 
solving the dynamical equation, an ISM density of n = 
1 cm is taken. For the initial explosion conditions, two 
sets of parameters are used. The first set is representative 
for GRBs: I take E = 10 52 erg and . M = 2 x l(r 5 M Q , 
resulting in T = 278 ([Huang et alJll999f ). The second 
set is representative for XRBs, and is chosen as follows. 
As an initial Lorentz fact or I take a fiducial value To = 3 
(|Miller- Jones et al.|[2005T) . Observed XRB radio blobs 
are typically emitted when the luminosity is close to the 
Eddington luminosity, and the flux rise time is a few 
days, hence the total energy release d is of the order of 
~ 10 45 erg (e.g.. lFender~etaTI l200l . With T = 3 this 
leads to an ejected mass of M — 3 x 10 23 gr. 

The results in Figures Q]- [3] are given for both the adi- 
abatic scenario and the radiative scenario. As expected, 
the results asymptote to the known solutions, which can 
be divided into 3 regimes. (I) Initially, r ~ Tq, and 
R oc t; (II) Tm ^> M, and r>l. In the adiabatic sce- 
nario, this leads to R(t) oc i 1/4 and T(t) oc t -3 / 8 , while 
in the radiative scenario, R(t) oc i 1 / 7 and T(t) oc i~ 3 / 7 . 
(Ill) For f) <C 1, in the adiabatic scenario R(t) oc £ 2 / 5 and 
j3(t) oc £ -3 / 5 , while in the radiative scenario, R(t) oc t 1 / 4 
and j3(t) oc t- 3 / 4 . 

4. SUMMARY AND DISCUSSION 

In this letter, I have revisited the dynamics of a blast 
wave propagating into an ISM, as is expected in many 
astronomical objects, such as GRBs, XRBs and Super- 
novae. I derived an equation ([5]) that determines the 
evolution of the blast wave Lorentz factor as the plasma 
collects material from the ISM and decelerates. Ana- 
lytical solutions in both the adiabatic and the radiative 
regimes are found in the asymptotic limits r » 1 and 
^ < 1 (i )2.ip . Numerical integration is easily carried, 
and the resulting dynamics valid in the full regime is 
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Fig. 2. — Temporal evolution of the momentum, T/3. All param- 
eters values are similar to Figure [l] 
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Fig. 3. — Temporal evolution of the shock radius, R. All param- 
eters values are similar to Figure [l] 

presented i]3l 

The dynamical equation is different than the 
dynamical equations derived earlier by several 



authors featz fc Piranl Il997t iHuang et all 119991: 
iChiang fc Dermerl 119991 : IPiranl 11999ft . This is due to a 
conceptual difference: earlier works considered a "toy 
model" , in which the basic assumption is that the 
interaction between the blast wave and the interstellar 
material (ISM) can be described by a series of inelastic 
collisions. As opposed to that, here I consider the full 
energy-momentum tensor, which takes into account 
the contribution of the collected material to both the 
energy and the pressure in the shocked region. Such an 
inclusion was neglected in earlier works. 

While the results of earlier works retrieve the correct 
asymptotic behavio r in the limits F 3> 1 and /3 <C 1 
(jHuang et al.l ll999). the evolution of the Lorentz factor 
derived in these works is different than the evolution de- 
rived here by a numerical factor of tens of percents. In 
the radiative scenario, the difference in the derived value 
of the Lorentz factor at a given radius is up to > 30%. 
This result is not surprising, given that the main dif- 
ference between the dynamics derived here and the dy- 
namics derived in former works lies in the inclusion of 
7 ~ 4/3 (for r ^> 1). Interestingly, also in the adiabatic 
scenario the numerical difference is larger than 12%. 

In recent years, a renewed interest in calculating the 
expected flux from the interaction of an expanding 
blast wave with its environment had emerged. Works 
had bee n carried in the context of GRB afte rglow 
emiss i on (IGranot fc Ramirez - Ruiz 20iot lGranot fc Piranl 
[Ml Ivan Eerten fc MacFadvenl I2012t IXu et all 12012). 
the evolution of the obs e rved radio blobs seen in XRBs 
(jShen fc Matznerl 120121: iNaravan fc McClintockJ I2012T) 
as well as the evolution of emission from supernovae 
(|Chakraborti fc Ravll20liT ). The dynamical calculations 
presented here are an obvious essential step in these cal- 
culations. 
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